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Abstract 



^ ■ We demonstrate that with a stepwise introduction of complexity to a model of an electron sys- 

tem embedded in a photonic cavity and a carefully controlled stepwise truncation of the ensuing 
many-body space it is possible to describe the time-dependent transport of electrons through the 
system with a non-Markovian generalized quantum master equation. We show how this approach 
retains effects of an external magnetic field and the geometry of an anisotropic electronic system. 
^2 ■ The Coulomb interaction between the electrons and the full electromagnetic coupling between the 

Q ■ electrons and the photons are treated in a non-perturbative way using "exact numerical diagonal- 

ly^ ! ization" . 
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I. INTRODUCTION 

Advances in techniques constructing and experimenting with quantum electro dynamic 
circuits have resulted in systems with very strong electron-photon coupling flHS]. Tradi- 
tionally, some version of the Jaynes-Cummings model [^ is used to describe the energy 
spectrum of the closed system or its time evolution j5| . Recently, we have shown that for 
a strong electron-photon coupling in a semiconductor nanostructure the Jaynes-Cumming 
model may not be adequate and one may have to consider a model with more than two elec- 
tron levels and the diamagnetic term in the coupling [6]. In continuation we have used our 
experience with describing the dynamics of open systems in terms of the generalized master 
equation (GME) [71, |8| to start the exploration of time-dependent transport properties of 
circuit quantum electrodynamic (circuit-QED) systems [9]. 

Here, we will describe our approach with a special emphasis on what we call: "A stepwise 
introduction of complexity to a model and a carefully controlled stepwise truncation of the 
ensuing many-body space" . We will discuss technical issues that are common to models of 
different phenomena and fields, but we will use the model of Coulomb interacting electrons 
in a photonic cavity as an example to display our approach and findings. 
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FIG. 1. A schematic view of the total system consisting of a finite quantum wire (green, the 
central electronic system) in a photon cavity (bold black), coupled to external left (red) and right 
(blue) leads at different chemical potentials, /iL and /iR. An external homogeneous magnetic field 
B (magenta) is perpendicular to the system and the leads. The lengths are not shown to scale. 
The length of the central system is 300 nm, but the characteristic lengths of the cavity are on 
the millimeter scale. The GaAs leads and central system contain a quasi-one-dimensional (QID) 
electron gas in 1 to 4 subbands of width 60 to 120 nm for the parameters to be introduced below. 
The cavity photon modes (yellow) are standing waves in the z-direction with an electric component 
in the x- or y-direction. The applied bias between the left and right leads is indicated by placing 
them at different levels with respect to the central system. 



II. THE MODEL OF THE CLOSED SYSTEM 



The closed electron-photon system is a finite quasi-one-dimensional (QID) quantum wire 
placed in the center of a rectangular photon cavity (see Figure [I]) . It is described by the 



Hamiltonian J9|] 

Ho = ^ Eid\di + hjJo)a + - y^^{ij\Vcou\\rs)d\d\dsdr 

i ijrs 

+ S,J24dJ9^J{a + a^} (1) 

where Ei is the single-electron spectrum for the finite quantum wire with hard walls at 
X = ±Lx/2 and parabolic confinement in the ^-direction with characteristic energy HQq. A 
static classical external magnetic field B = V x Agxt = Bz renormali zes the freq uency of 
the ^/-confinement Q^^ = uj'^ + Q^ and the natural length scale a^ = -\//z/(m*^^) with the 
cyclotron frequency uJc = eB /{rrfc). di is an annihilation operator of the non-interacting 
single-electron state (SES) \i) with energy Ei^ and a is the annihilation operator for the 
single-photon mode with energy hw. The kernel for the Coulomb interaction of the electrons 
is 

with the small regularization factor selected such that ry/a^ ?^ 7.1 x 10~^ when a^ '^ 33.5 
nm at 5 = 0.1 T for GaAs parameters, i. e. the eflFective mass m* = O.OGTmg, and the 
dielectric constant k = 12.4. The second line in the Hamiltonian ([T]) is the paramagnetic 
interaction between electrons and photons (— J rfr j • A)/c, and the last line stems from the 
diamagnetic term in the interaction (— e Jrfr pA^)/{2rrfc^). In terms of the field operators 
the charge current density and charge density are, respectively, 

{^t (tt^) + (7r>t) ^} ^ (3) 



27TZ* 

and 

p = -67/;^, (4) 

where 

^ = (p + - Aext) . (5) 

The photon cavity is assumed to be a rectangular box (x,y,z) G {[— ac/2, ac/2] x 
[— ac/2,ac/2] x [— rfc/2, (ic/2]} with the finite quantum wire centered in the z = plane. 
In the Coulomb gauge the polarization of the electric field can be chosen parallel to the 
transport in the x-direction by selecting the TEqh mode, or perpendicular to it by selecting 
the TEioi mode 

A is the amplitude of the cavity vector field defining a characteristic energy scale £c = 
eA^w^w/c = g^^ for the electron-photon interaction and leaving an eflFective dimensionless 
coupling tensor 

9^J = '^jdr [^*(r) {(e • tt) ^,-(r)} + {(e ■ tt) ^,(r)}* ^.-(r)], (7) 



defining tlie coupling of individual single-electron states \i) and \j) to the photonic mode. 
In the calculations of the energy spectrum of the Hamiltonian ([T]) we will retain all resonant 
and antiresonant terms in the photon creation and annihilation operators so we will not use 
the rotating wave approximation, but in the calculations of the electron-photon coupling 
tensor ([7j) we assume a^, Lx « a^ and approximate cos(7r{x, y}/ac) ~ 1 in Eq. (|6j) for the 
cavity vector field A. 

The energy spectrum and the states of the Hamiltonian for the closed electron-photon 
system have to be sought for an unspecified number of electrons as we want to open the 
system up for electrons from the leads later. We will be investigating systems with few elec- 
trons present in the finite quantum wire, but it is not a trivial task to construct an adequate 
many-body (MB) basis for the diagonalization of Hq since in addition to geometrical and 
bias (set by the leads) considerations we have strong requirements set both by the Coulomb 
and the photon interaction. Our solution to this dilemma and a mean to keep a tight lid 
on the exponential growth of the size of the required many-body Fock space is to do the 
diagonalization in two steps. 

First, we select the lowest A^ses single-electron states (SESs) of the finite quantum wire. 
These have been found by diagonalizing the Hamiltonian operator for a single electron in 
the QID confinement and in a perpendicular constant magnetic field in a large basis of 
oscillator-like wave functions. Originally, we constructed a many-electron Fock space with 
^MES = 2^sEs states |/x) [lO|, lH]. (We use Latin indices for the single-electron states and 
Greek ones for the many-electron states). This "simple binary" construction for the Fock- 
space does not deliver the optimal ratio of single-, two-, and higher number-of-electrons 
states for an interacting system when their energy is compared. Usually one ends up with 
too few SESs compared to the MESs. Here we will select 18 SESs and construct all possible 
combinations of 2-4 electron states. This can be refined further. These MESs, which are in 
fact Slater determinants, and which we denote as \fi) (with an angular right bracket) are then 
used to diagonalize the part of the Hamiltonian ([T]) for the Coulomb interacting electrons 
only, supplying their spectrum E"^ and states |/i) (denoted now with a rounded right bracket). 
The eigenvectors from the diagonalization procedure define the unitary transform between 
the two sets of MESs, |/x) = X^^^ V^a|a). As the action of the creation and annihilation 
operators is only known in the non-interacting electron basis {|/i)} we need this transform 
to write Hq in the new Coulomb interacting basis {|/x)} 

Ho = Yl Im)^m(mI + f^^^^ + ^c 5Z \l^)W^4djV\iy){iy\ g^J [a + o)] 

+ ^^ i-w) ^ \^'W^^4d^VW){u\ I ba + 1) + 1 (aa + a^a^) \ . (8) 

In order to finally obtain the energy spectrum of the electron-photon Hamiltonian ([8]) we 
need to construct a MB-space |/x) ® |Aph) — > |M)e-EM out of the Coulomb interacting 
MESs l/i) and the eigenstates of the photon number operator |Aph). To properly take 
account of the effects of the Coulomb interaction we selected a large basis {|/i)} and many 
photon states since the strong coupling to the cavity photons requires many states. For 
the system parameters to be introduced later we find that basis build up of the 64 lowest 
in energy Coulomb interacting MESs and 27 photon states is adequate for the transport 
bias windows and the electron-photon coupling to be selected later. The catch is that the 
unitary transform has to be performed with the full untruncated basis since V can not be 



truncated. In other words each of the 64 interacting MESs |/x) remains a superposition 
of all of the 2^^^^ noninteracting MESs |//) with the same number of electrons N^^ i. e. 
l/i) is a linear combination of a subset of ( ^^^) terms. (A similar issue is met when the 
dimensionless coupling tensor gij ([71) is transformed from the original single-electron basis 
to the single-electron states \i)). 

The diagonalization of Hq produces the new interacting electron-photon states \jll) = 
Xla ^i2a\^)e-EM with a kuowu integer electron content, but an indefinite number of photons, 
since the photon number operator does not commute with Hq. The MB energy spectra for 
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FIG. 2. The single-electron energy spectrum of the leads (left panel) versus the "subband mo- 
mentum" qa^, and the many-body energy spectra for Coulomb interacting electrons coupled to 
quantized cavity photon modes with the electric component polarized along the finite quantum wire 
(x-polarization, center panel), and perpendicular to the wire (y-polarization, right panel) versus 
the electron-photon coupling strength g^^ — £c- In the energy range shown there are states with 
no electrons (horizontal solid green), one electron (solid red), and two electrons (dashed blue). In 
order to reach convergence for high S^ in this figure we use A/'ses = 200 here. S = 0.1 T, hTt^) = 1.0 
meV, hijj = 0.4 meV, W}^ = 1.0 meV, L^ = 300 nm, m* = O.OGTmg, /^ = 12.4, the dielectric 
constant of GaAs. 



the X- and y-polarization of the photon mode are shown in Fig. [2] in comparison with the 
subbandstructure of the single-electron spectrum of the leads. The diflFerence in the energy 
spectra for the x- and the y-polarization stems from the anisotropic electron system. The 
photon mode selected with energy 0.4 meV is close to be in resonance with the motion in 



the x-direction, but quite far from the fundamental energy in the ^-direction, 1.0 meV. The 
spectrum for the x-polarization thus displays stronger dispersion and interaction of levels 
than the spectrum for the ^-polarization. 

The states of the closed system in Fig. [2] have a definite number of electrons and an 
undetermined number of photons of one given polarization. Later when the system is coupled 
to the leads one can expect the electrons entering the system to radiate photons with any 
polarization. A preferred polarization could though be infiuenced by the aspect ratio of the 
cavity. In the treatment of the time-dependent system to follow we assume the electrons 
only to radiate photons with one given polarization. This is done to facilitate the numerical 
calculations and the analysis of a system with a strong spatial anisotropy. 



III. OPENING OF THE SYSTEM, COUPLING TO EXTERNAL LEADS 

At time t = the closed system of Coulomb interacting electrons coupled to cavity 
photons described by the Hamiltonian (|8j) is connected to two external semi-infinite quasi- 
one-dimensional quantum wires in a perpendicular magnetic field J8|. We use an approach 
introduced by Nakajima and Zwanzig to project the time evolution of the total system onto 
the central system by partial tracing operations with respect to the operators of the leads 
(12L ll3|. The coupling Hamiltonian of the central system (i. e. the short QID wire) to the 
leads is of the form 

Hrit) = J2x'{t) fdq {Ti4,d, + (Ti)*4c,,} , (9) 

i,l -^ 

where / G {L, R} refers to the left or the right lead, and xK't) is the time-dependent switching 
function of the coupling. The operators Cqi and c\ annihilate and create an electron in the 
/-lead with a quantum number q referring both to the continuous momentum q and the 
subband n\ see Fig. [2] for the corresponding energy spectrum. To represent the geometry 
of the leads and the central system, the coupling tensor T^- of single-electron states \q) in 
the lead / to single-electron states states \i) in the system is modeled as a non-local overlap 
integral of the corresponding wave functions in the contact regions of the system, f^^^, and 
the lead /, f^/[Il] 

T!,= f rfrrfr^K(rO]>f(r)^^r,rO. (10) 



r^x^^Li 



The function 



i(r,r') = ^^exp [-S[ix - x'f - 5[{y - y'f] x exp (^^^^^) (H) 

with r G f^g and r' G ^i defines the 'nonlocal Gaussian overlap' determined by the constants 
5i and ^2, and the affinity of the states in energy A^. The energy spectra of the leads are 
represented by e\q). 

The time-evolution of the total system is determined by the Liouville-von Neumann 
equation 

ihW(t) = [Hit), W{t)i W{t < 0) = plPrPs, (12) 



with W the statistical operator of the total system and pi the equilibrium density operator 
of the disconnected lead / G {L^ R} having chemical potential pi 

^^ ^ Tr^{e-^(^^-/^^^^)}' ^^^^ 

where Hi is the Hamiltonian of the electrons in lead / G {L, R} and Ni is their number 
operator. The Liouville-von Neumann equation f IT2l) is projected on the central system of 
coupled electrons and photons by a partial tracing operation with respect to the operators 
of the leads. Defining the reduced density operator (RDO) of the central system 

ps(t) = TrLTrRH^(t), ps(0) = ps, (14) 

we obtain an integro-differential equation for the RDO, the generalized master equation 

(GME) 



Ps{t) = -'-[Hs^Ps{t)]-^ 



HT{t), f M [U{t - t')Hr^{t')U^ {t - t'), 
Jo 

Uo{t - t^)ps{t^W{t - t>LPR] ] } , (15) 

where the time evolution operator for the closed systems of Coulomb interacting electrons 
coupled to photons on one hand, and on the other hand noninteracting electrons in the leads 
is given by U{t) = exp {—i^H^ + i^coui + ^em + ^l + H^)t/h}^ without the coupling to the 
leads H^{t). Here, Hq is the Hamiltonian for the electrons in the central system, i^coui 
their mutual Coulomb interaction, i^EM is the Hamiltonian for the photons in the cavity 
together with their interaction to the electrons, and i/L,R are the Hamiltonian operators for 
the electrons in the left and right leads. The time evolution of the closed system of Coulomb 
interacting electrons interacting with the photons is governed by U{){t) = exp {—iHot/h}. 
The GME (IT5l) is valid in the weak system-leads coupling limit since we have only retained 
terms of second order in the coupling Hamiltonian H^ in its integral kernel. It should though 
be stressed that we are not approximating the GME to second order in the coupling as its 
integral structure effectively provides terms of any order, but with a structure refiecting the 
type of the integral kernel. 

Commonly, the GME is written in terms of spectral densities for the states in the system 
instead of the coupling tensor (ITOl) . We do not make this transformation but the spectral 
densities for the 10 lowest SESs used in our system have been presented elsewhere [9] (Here 
the overall coupling is one quarter of the value used in Ref. [9]). The spectral density is a 
particularly useful physical concept as it demonstrates the spectral broadening of the SESs 
in the finite system due to the coupling to the leads. 

IV. TRANSPORT CHARACTERISTICS 

The time-dependent coupling betweeen the leads and the central system is modeled by 
the switching functions xK't)- These functions may be considered input elements of the 
transport problem: stepwise functions, periodic, relatively phase-shifted, etc. In the follow- 
ing examples the left and right leads are coupled simultaneously smoothly to the central 
system by use of the switching function 



x'W=(l-^^)' 1^{L,R} 



(16) 



with Qf^ = 0.3 ps ^. The temperature of the leads T = 0.5 K, and the overaU couphng 
strength g^aj = 13.3 meV, is much lower than in our earlier calculation [9]. We choose 
6[ 2a^ = 0.4916 meaning that states of a lead and the central system with considerable 
charge density within a length equivalent to 2a^ (a^ ^ 33.5 nm here) could be well coupled. 
The energy of the photon mode is huj = 0.4 meV, and the electron confinement in 
the ^-direction has the energy scale HQq = 1.0 meV. The energy separation of the lowest 
states for the motion in the x-direction is lower, ^ 0.2 meV. The Coulomb interaction has a 
characteristic energy scale 0.5e^/(/^a^) ^ 1.7 meV. The external magnetic field is low enough 
so that only the highest lying SESs show any effects of the Lorentz force. We expect thus 
the transport properties of the system to be anisotropic at a low energy scale with respect 
to polarization of the photon field along (x-direction) or perpendicular (^-direction) to the 
transport. 
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FIG. 3. The total mean number of electrons (A^e(^)) (left panel), and the total mean photon 
number (A^ph(^)) (right panels) as a function of time, polarization of the photon field (x or y) and 
bias window low (1: /j^l = 2.0 and jir = 1.4 meV) or high (h: /j^l = 3.0 and jir = 2.5 meV). B = 0.1 
T, ^™ = 0.1 meV, hQo = 1-0 meV, L^ = 300 nm, huj = 0.4 meV, A^ = 0.25 meV, gl^a^^ = 13.3 
meV, S[2al = 0.4916, m* = 0.067me, and k. = 12.4. 

In the present calculations we start with a central system empty of electrons, but with 
one or four photons in the cavity. For A^ses = 18 the state with no electron, but one photon, 
is 1 2) for both polarizations. The state containing 4 photons and no electrons is |14) for 
the x-polarization and |15) in the case of the ^-polarization. The time evolution of the 
total mean electron number (A^e(^)) is seen in the left panel of Fig. [3] and the total mean 
number of photons {N^h{t)) is displayed in the lower right panel for the case of initially 
one photon and in the upper panel for initially 4 photons. The charging of the system is 
in most cases fastest for the higher bias window as the active states of the central system 
are well coupled to states in the leads that carry a large current in this energy range. The 
presence of photons in the system sharply diminishes the charging speed, especially for their 
polarization along the transport direction (x-direction). The time evolution of the total 



8 



mean number of photons in the system does not give much insight into what is happening 
in the system, but it is though clear that it varies more in case of the x-polarization. 

Very similar information can be read from the graphs of the total currents in the left 
(L) or the right (R) leads shown in Fig. [H The negative values for the current into the 
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FIG. 4. The total current from the left lead (L), and into the right lead (R) as a function of time 
for g^^ = 0.1 meV and polarization of the photon field (x or y). Initially, at t = there is one 
photon in the cavity (top panel), or 4 photons (bottom panel). The bias window is low (/x^ = 2.0 
and fiR = 1.4 meV) for the left panels and high (/Xi, = 3.0 and /j.r = 2.5 meV) for the right 
panels. B = 0.1 T, h^o = 1.0 meV, L^ = 300 nm, hw = 0.4 meV, jil = 2.0 meV, jir = 1.4 meV, 
A^ = 0.25 meV, gl^a^'^ = 13.3 meV, 5^ 2^1 = 0.4916, m* = 0.067me, and /^ = 12.4. 

right lead indicate a flow from the lead to the central system. A note of caution here is 
that we are concentrating on the charging time-regime here and we are not trying to reach 
a possible steady state, but clearly we are most likely in a Coulomb blocking range. It is 
our experience that a considerable probability for two electrons in the system is only seen 
after several two-electron states are in or under the bias window. The reason for this is most 
likely the different coupling between states in the leads and the system, how far the system 
is from equilibrium and, and how the coupling strength influences the rate of occupation 
of various dynamically correlated states. We have verifled that a still higher bias leads to 
a nonvanishing steady state current. Here, we always have at least one photon initially in 
the cavity and we notice that the charging rate and the currents are mostly higher for the 
y-polarization. The time-scale for the charging in the x-polarization gets very long as the 
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photon number is increased from 1 to 4. This is not a total surprise since the energy of the 
photons is closer to characteristic excitation energies for a motion in the x-direction. 

A detailed view of the charging processes can be obtained by observing the time- 
dependent probabilities for occupation of the available many-body states (MBS) by electrons 
or photons. In Fig. [5] we see the mean charge in the MBS for the transport through the 
lower bias window jlll = 2.0 and /llr = 1.4 meV. The MBS \JLt) are numbered according to 
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FIG. 5. The mean number of electrons (A^e(^)) in a MBS \fl) for a low bias window {/j^l = 2.0 and 
/ii^ = 1.4 meV) for x-polarization (left) and y-polarization (right) as a function of time. Initially, 
at t = 0, there is one photon in the cavity (top panels), or 4 photons (bottom panels). The initial 
number of electrons is zero in all cases. S = 0.1 T, g^^ = 0.10 meV, hcj = 0.4 meV, hfto = 1.0 
meV, A^ = 0.25 meV, gl^a^J^ = 13.3 meV, S^^al = 0.4916, L^ = 300 nm, 
A^ = 12.4. 
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increasing energy. We see clearly that the system is not close to equilibrium and the charge 
is "scattered" to more states for the x-polarization than the ^-polarization. The presence 
of photons has large effects on the electrons in the system. 

Very similar story can be said about the results for the higher bias window, //^ = 3.0 
and fiR = 2.5 meV displayed in Fig. [61 Especially interesting is to see that the results for 
both polarizations are almost identical for the higher and lower bias window for the case of 4 
photons initially in the central system. The effects of the bias window are washed out by the 
strong interaction with the quantized electromagnetic field in the cavity. We have though 
to bring in a note of caution here and admit that this effect should be studied further using 
a larger MB basis for the RDO. 

A complete picture can not be reached without exploring the time-evolution of the mean 
photon number per MBS presented in Fig. [7] for the lower bias window and in [8] for the 
higher bias window. Again we notice that the state of the system, the distribution of the 
photon component into various MBS, is almost independent of the the bias window in the 
case of 4 photons initially in the cavity. 
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FIG. 6. The mean number of electrons (A^e(^)) in a MBS \fl) for a high bias window (jil — 3.0 and 
fiR = 2.5 meV) for x-polarization (left) and ^/-polarization (right) as a function of time. Initially, 
at t = 0, there is one photon in the cavity (top panels), or 4 photons (bottom panels). The initial 
number of electrons is zero in all cases. S = 0.1 T, g^^ = 0.10 meV, huj = 0.4 meV, h^o = 1.0 
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But, here another very important fact about the system evolution becomes evident. The 
occupation of the initial photon state seems to vary much faster with time than the total 
mean number of photons shown in Fig. [31 The slow decay of the charging current in Fig. 
m for the system with 4 photons initially might also indicate that radiation processes here 
are slow. Why then do we have a fast redistribution of the photon component between 
the available MBS initially, during the switch-on process? The resolution of this dilemma 
comes from remembering the structure of the interaction terms. Part of the interaction is an 
integral over the term j • A. We have a vector field initially present in the system. The high 
current into the system during the transient phase, when the contacts between the leads 
and the central system are switched on, creates a strong interaction that can "scatter" the 
electrons and the photons to different MBS, few or many, depending on selections rules and 
resonances. This initial rushing of electrons from both ends of the finite quantum wire is 
a longitudinal (irrotational) current that enhances the coupling to the photon field. This 
eflFect can be seen very well by comparing the photon distribution into the MBS at the low 
bias window and for the case of only one photon initially present in the cavity. 

In case of the x-polarization we see that the system seems to develop a higher "impedance" 
acting against its charging as the number of photons is increased, and at the same time the 
electrons that enter the system are scattered to a wide range of MBS. The "impedance" 
has to be understood from the fact that the photon field is polarizing the electron density 
inside the system along the x-direction. Effectively, the states may gain coupling strength 
to the leads, but that works in both ways. The electrons enter the system easily and leave it 
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FIG. 7. The mean number of photons (A^ph(^)) in a MBS \fl) for a low bias window (fi^ = 2.0 and 
fiR = 1.4 meV) for x-polarization (left) and ^/-polarization (right) as a function of time. Initially, 
at t = 0, there is one photon in the cavity (top panels), or 4 photons (bottom panels). 5 = 0.1 
T, ^™ = 0.1 meV, hw = 0.4 meV, hfto = 1-0 meV, A^ = 0.25 meV, g^aH^ = 13.3 meV, 
S{^^al = 0.4916, L^ = 300 nm, m* = 0.067me, and n = 12.4. 



equally easily at the same end before they can be transported through it. The charging of the 
system is counteracted by an effective scattering of a localized photon field. A phenomena 
that should be similar for a system with localized phonons. 

V. SUMMARY 

In this publication we have shown how we have successfully been able to implement a 
scheme which we wish to call: "A stepwise introduction of complexity to a model description 
and a careful counteracting stepwise truncation of the ensuing many-body space" to describe 
a time-dependent transport of Coulomb interacting electrons through a photon cavity. We 
have used this approach to demonstrate how the geometry of a particular system leaves its 
fingerprints on its transport properties. We guarantee geometrical dependence by using a 
phenomenological description of the system-lead coupling based on a nonlocal overlap of the 
single-electron states in the leads and the system in the contact area, and by using a large 
basis of SESs in order to build our MBS. The Coulomb interaction is implemented by an 
"exact numerical diagonalization" , and the coupling to the quantized electromagnetic cavity 
field of a single frequency is carried out using both the paramagnetic and the diamagnetic 
part of the charge current density. The coupling of the electrons and the photons is also 
treated by an exact numerical diagonalization without resorting to the rotating wave ap- 
proximation. This approach is thus applicable for the modeling of circuit-QED elements in 
the strong coupling limit. 
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FIG. 8. The mean number of photons (A^ph(^)) in a MBS \fl) for a high bias window (fiL = 3.0 and 
fiR = 2.5 meV) for x-polarization (left) and ^/-polarization (right) as a function of time. Initially, 
at t = 0, there is one photon in the cavity (top panels), or 4 photons (bottom panels). 5 = 0.1 



T, 

6{ -2 



^EM 



az, = 0.4916, L^ 



0.1 meV, hw = 0.4 meV, h^o = 
300 nm, tti* = O.OGTTTig 



1.0 meV, A^ = 0.25 meV, gl^ati = 13.3 meV 



and hi = 12.4. 



We have deployed the GME method to include memory effects without a Markov approx- 
imation. Moreover, the GME approach is utilized to describe the coupling of the electron- 
photon system to the external leads. At this moment in the development of the model we 
ignore the spin variable of freedom, but it can easily be included. The inclusion of the spin 
is (more or less) straight forward, but it requires a double amount of computer memory and 
an increased number of MBSs in the GME solver. This is work in progress now. 

We concentrate our investigations on the charging regime of the system and find that it 
is strongly infiuenced by the presence of photons in the cavity, their polarization and the 
geometry of the system. We consider these results presented here as the mere initial steps 
into the exploration of a fascinating regime of circuit-QED elements. 

For the numerical implementation we rely heavily on parallel processing, but we foresee 
further refinement in the truncation schemes for the many-body spaces and in the paral- 
lelization that will allow us to describe systems of increased complexity. 
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